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The Mott insulating state formed from bosons is ubiquitous in solid *He, cold atom systems, 
Josephson junction networks and perhaps underdoped high-Tc superconductors. We predict that 
close to the quantum phase transition to the superconducting state the Mott insulator is not at all 
as featureless as is commonly believed. In three dimensions there is a phase transition to a low tem- 
perature state where, under influence of an external current, a superconducting state consisting of a 
regular array of 'wires' that each carry a quantized flux of supercurrent is realized. This prediction 
of the "type-II Mott insulator" follows from a field theoretical weak-strong duality, showing that 
this 'current lattice' is the dual of the famous Abrikosov lattice of magnetic fiuxes in normal super- 
conductors. We argue that this can be exploited to investigate experimentally whether preformed 
Cooper pairs exist in high-Tc superconductors. 
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I. INTRODUCTION 

The Yin- Yang mystique in Asian philosophy has found 
a remarkably literal incarnation in modern physics in the 
form of the duality principle^. An elementary example 
of the idea that 'opposites form a unity' is the particle- 
wave duality of quantum mechanics. This was surpassed 
by the identification of the Kramers-Wannier or weak- 
strong duality structures in quantum field theory, even- 
tually leading to the rich dualities of string theorjEl. A 
pedestrian example that will play a role in the back- 
ground of the present story is the electromagnetic du- 
ality, stating that in a world where magnetic monopoles 
have a similar standing as electrical charges, the 'oppo- 
site' electrical and magnetic universes are mathematical 
mirror imaged. 

There is yet a deeper level that becomes particularly 
explicit when dealing with the strongly-interacting quan- 
tum many-particle systems of condensed matter physics. 
Such systems will typically have ordered ground states 
breaking symmetry spontaneously, for instance the su- 
perconducting state. With mathematical topology it is 
then possible to identify field configurations that are 
uniquely associated with the restoration of the bro- 
ken symmetry: the topological excitations, such as the 
Abrikosov vortices in a superconductor. Upon increasing 
the quantum fluctuations of the collective state (e.g. in- 
creasing charging energy in the superconductor) at some 
point the system will undergo a zero-temperature quan- 
tum phase transition (QPT) where the system 'melts' 
into a quantum disordered state*^. The weak-strong 
duality principle now prescribes that this quantum dis- 
ordered state can always be viewed as some ordered state 
formed from the topological excitations associated with 
destroying the order of the ordered state. 

The archetypical example is the "vortex-boson" or 
"Abelian-Higgs" duality in two space and one time di- 
mensions (2-l-lD), associated with a system of interacting 
bosons living on a lattice, undergoing a superconductor- 



insulator QPT. The simplest microscopic model of rele- 
vance to this physics is the Bose-Hubbard HamiltoniarP^l, 
realized literally in cold atom systems for t he neutral 
casC"'^ and Josephson junction networkJi^l^ for electri- 
cally charged bosons. For neutral bosons it reads as 

i?BH = -t'Y.'y^l^j + + U^ni{ni - 1) - ni, 

{ij) i i 

(1) 

which can be straightforwardly extended to the charged 
case by coupling in the electromagnetic gauge fields. We 
specialize to the case with an integer number of bosons 
per lattice site ("zero chemical potential"). For small 
charging energy U the bosons will condense into a super- 
fluid/superconductor. However, when U/t w 1 a quan- 
tum phase transition occurs to a Bose-Mott insulator. 
The charging energy exceeds the kinetic energy with the 
effect that the bosons localize: a Mott gap opens and the 
lowest lying excitations are the doublons (extra boson) 
and holons (missing boson). This Bose-Mott insulator is 
conventionally considered to be a completely featureless 
state, not breaking any symmetry. 

However, the naive picture of localized bosons is 
fiawecP, since the ground state energy will always be 
lowered by virtual exchange processes. Moreover, al- 
though often not realized, the true nature of the Mott in- 
sulator clo se to the QPT is revealed by the vortex duality 
perspectiv d^ l ^^* ^^. In 2-l-lD vortices are point particles, 
and they embody the virtual quantum fluctuations in 
the superconductor as closed loops of vortex-antivortex 
worldlines. These loops grow in size when approaching 
the QPT, to 'blow out' at the transition and the Mott 
insulator corresponds with a tangle of free vortex and 
antivortex worldlines (Fig. [2]). Elegantly, the vortex- 
vortex interactions can be parametrized in terms of ef- 
fective U{1) gauge fields, and this tangle of worldlines is 
therefore identical to a relativistic (Higgs) superconduc- 
tor, where the Higgs mass is coincident with the Mott 
gap, while the holon and doublon excitations have the 
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same status as the 'massive photons' of this dual vor- 
tex superconductor^. This suggests that there is more 
going on than the featureless state one infers from the 
strongly-coupled, atomic limit 'canonical view'. 

However, one has to now consider the thermodynamics 
of the dual superconductor. The vortex superconductor 
is charged and therefore the interactions between the dual 
vortices in the dual condensate are short-ranged. Since 
these are particles in 2+ ID, they will proliferate at any 
finite temperature. There is therefore no thermodynam- 
ical phase associated with the dual superconductor: one 
recovers the featureless Mott insulator. This situation is 
drastically different in three space dimensions. 

Until very recently it was not quite known how to for- 
mulate vortex duality in the natural 3+1 dimensions of 
the physical world. The obstruction was of a techni- 
cal nature. In three space dimensions vortices are lines, 
which implies that in spacetime (as of relevance to the 
zero-temperature physics) vortices correspond to quan- 
tum strings (Fig. [T]). Instead of the tangle of vortex 
worldlines in 2+lD, the dual condensate now consists of 
a 'foam' formed from the vortex worldsheets. Although 
unrelated to fundamental string theory, one cannot rely 
on the standard methods of qua ntum field theory for the 
description of such a condensat^SEII (ge e also Ref. [T5)) . 
Its workings were tackled only recentljff^. We will review 
this below, but the outcome is actually rather straight- 
forward: this 'two-form Higgs phase' is qualitatively very 
similar to the standard (relativistic) superconductor, the 
main difference being in the counting of degrees of free- 
dom. 

Here we report on the extension from the neutral su- 
perfluid to the charged superconductor in 3-l-lD. Al- 
though far from self-dual, the charged Bose-Mott insu- 
lator as a 'dual stringy superconductor' is behaving as a 
normal superconductor to the extent that the topology 
of the phase diagram of the Mott insulator in 3+lD is 
a 'dual mirror image' of the phase diagram of a normal 
superconductor: our main result. Fig. [4] In a normal 
superconductor, the control parameters are temperature 
and applied magnetic field. As we will discuss in de- 
tail, the magnetic field dualizes into applied current in 
the Mott insulator, and after re-identification of this axis 
the phase diagrams on both sides of the superconductor- 
Mott insulator transition acquire the same topology. To 
read off the physics of the Mott insulator, one can just 
depart from the standard wisdoms for superconductors 
using the 'dual dictionary' summarized in Table [l] 

A quite counterintuitive prediction follows: upon re- 
ducing temperature, one will find a thermal phase tran- 
sition to the Mott insulator with the same thermody- 
namical {XY) signatures as in a bosonic (local pair) su- 
perconductor. In a normal superconductor one applies 
magnetic fields to probe the "generalized rigidity"^ of 
the ordered state. The dual of the magnetic field be- 
comes in the Mott insulator the electric current. Just 
as in the Meissner phase the magnetic field is expelled, 
in the "type-I Mott insulator" the electric current is ex- 
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photon field 
applied magnetic field 
London penetration depth 
flux quantum 
Meissner state 
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type-II Mott insulator 

vortex condensate |<1?|'^ 

dual gauge field b 

applied current J 

Mott proximity depth Am 

current quantum Iq 
insulating state 
current line lattice 
superconductor 
Maxwell vacuum 



TABLE I: Duality dictionary. The superconductor and the 
Mott insulator viewed as vortex superconductor are quite 
similar, but the meaning of the various physical quantities 
'turns upside down'. Magnetic fields turn into currents with 
the ramification that the Abrikosov lattice of magnetic fluxes 
penetrating the superconductor turns into a lattice of current 
fluxes penetrating the Mott insulator, the type-II phase. The 
'mirror image' is not perfect (not self dual), which leaves more 
states to play with in the case of the Mott insulator than for 
the normal superconductor (see Sec. |VI[|. 



pelled (with an associated dual penetration depth): this 
is just showing that the system is an insulator. 

However, starting from local pairs with a very short 
coherence length, one is generically dealing with type-II 
behavior, both in the normal and dual superconductors. 
In duality language, the difference between type-I and 
type-II behavior is due to whether the disordering parti- 
cles/strings (vortices) have net attractive or repulsive in- 
teractions. In a normal superconductor a magnetic field 
that exceeds the lower critical field will penetrate in the 
form of an Abrikosov lattice of vortex lines, carrying each 
a quantized magnetic flux. The dual of the type-II super- 
conductor is the "type-II Bose-Mott insulator"!^ when 
the external current exceeds a "lower critical current" it 
will penetrate the Mott insulator in the form of a lattice 
of 'wires' carrying each a quantized supercurrent! Macro- 
scopically it will behave just as a normal superconductor, 
which turns into a dissipative metallic state at the ther- 
mal transition where the dual order disappears. To find 
out whether such a superconductor is actually a Mott in- 
sulator in disguise one has to design experiments which 
are the 'current analogues' of the decoration experiments 
that led to the discovery of the Abrikosov lattice. 

As an immediate application of the idea we suggest 
to search for type-II Mott insulation behavior in under- 
doped cuprate superconductors^widely believed to be 
dominated by phase fluctuations^^J. Many researchers in 
the field are by now convinced that the so-called pseudo- 
gap regime consists of pre-formed Cooper pairs (bosons) 
that bind at a higher temperature T* , whereas phase co- 
herence and hence superconductivity set in only at at 
lower temperature T^.. In this scenario, the transition 
from the superconducting to pseudogap phase is precisely 
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(a)vortex worldline (b)vortex worldsheet 



FIG. 1: Vortex excitations in spacetime. (a) A vortex particle in two spatial dimensions traces out a worldline in spacetime. It 
is parametrized by the line element J^{x). |(b)| A vortex line in three spatial dimensions traces out a worldsheet in spacetime. 
Now we need two indices to define a surface element of the worldsheet, as such it is parametrized by the two-form field J^x- 
The vortices interact by exchanging gauge fields, that couple locally to the worldline or worldsheet (b^x). 



of the Xy-disordering type handled so well by vortex du- 
ality (see also Ref . I^T|) . Therefore we propose that in the 
vincinity of the quantum phase transition, the pseudo- 
gap phase is a type-II Bose-Mott insulator, and suggest 
several experimental setup that may verify the formation 
of quantized current lines in Fig. [5j 

The paper is organized as follows. We briefly recollect 
how the Bose-Hubbard model at zero chemical poten- 
tial maps to the XY-mode\ in Sec. |ll] Also the well- 
established vortex dualization procedure for 2+1 dimen- 
sions is reviewed in Sec |III| Subsequently we will sum- 
marize the results of Ref. [Klin Sec. |IV[ where we derive 
the vortex duality in 3-l-lD as well, and show how to 
incorporate the electromagnetic field to model the su- 
perconductor in which we are interested. The remainder 
contains new results. From the vortex duality in 3-l-lD, 
the prediction of a dual Meissner effect for current and in 
particular quantized lines of electric current immediately 
follows in Sec. |V] The main result is the phase diagram 
in Sec. |Vlj We propose several experiments that may 
confirm the existence of the quantized current lines in 
Sec. IVIII . The conclusions in Sec. IVIIII are followed 
by Appendix |X] which discusses the possible relevance of 
our findings to the so-called "giant proximity effect" , and 
Appendix |B] regarding Lorentz transformations of vortex 
worldsheets. 



II. BOSE-HUBBARD MODEL 

Here we recall shortly the Bose-Hubbard model of 
bosons hopping on a hypercubic lattice in D dimensions. 
For more detailed work see Refs. I4|10|I5I Let us first de- 
rive the relativistic continuum Ginzburg-Landau model 
of a superconductor by straightforward coarse graining. 
The Hamiltonian of the Bose-Hubbard model is, 

(ij) i i 

(2) 



Here b\ and are boson creation and annihilation oper- 
ators that satisfy the commutation relation [b^ , fej] — (5,j-. 

The number operator is rii = b\b^. Furthermore, the en- 
ergy scales are the boson hopping t, the on-site repulsion 
U and the chemical potential /j,. We shall assume that 
the chemical potential is tuned so that there is an integer 
number no of bosons per site ( "zero chemical potential" ) . 
Then we can make a change of variables b\ = y/n^e^'^' , 
so that the new conjugate variables satisfy the commuta- 
tion relation [93^,71^] — iSij. Substituting this definition 
in Eq. ^ leads to, 

H = -J ^(1 - cos(^i - Lpj)) + ^ ^(ni - l)n^. (3) 

{ij) i 

Here we have defined J = tno and added an unimportant 
constant. The physics of the weak and strong coupling 
limits is immediately clear: for large J/U , we have a 
superfluid where spatial fiuctuations in the phase are 
very costly. For small J /U the on-site repulsion domi- 
nates and the bosons are confined to their lattice sites: 
the Mott insulator. 

For the quantum field-theoretic formulation, we move 
from a Hamiltonian to a Lagrangian formalism, by noting 
that the canonical momentum is tTj = huj, which leads 
to the Lagrangian by Legendre transformation (where 

a ,„ OH U ^ \ 
i 

i {i-j) 

Now we can take the continuum limit in D space di- 
mensions I— J d^x; Lpi — ifj — aViy9(a;), where a 
is the lattice constant. This leads to the partition tunc- 
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tion Z = e R'^'^ in imaginary time t = ir where, 

^ I drd^x \ja^-^ [-^{dr^Y - (V^)^]. (5) 

This is to be compared with the quantum action for a 
superfluid [cf. Eq. (3.13) in Rcf. Hn], 

5e = / drd^x [ - \h^n{drvf - \ti'P^{Vvf] . (6) 
J '-2 2 m 

Hence we identify the compressibility n = jj^, the su- 
perfluid density divided by the boson mass ^ = ^^^rj — 
and the superfluid velocity Cph = | VC/ J. The energy 
scale VUJ will play an important role in the discussion 
of the quantum of electric current later on. Defining the 
covariant derivative 9^'' = {-^dr,'^), we find a conve- 
nient form of the action, 



drd^x 



-Ja 



2-D 



(7) 



We are interested in charged superfluids, i.e. super- 
conductors where the bosons must couple minimally to 
the electromagnetic potential, or photon field. Recall 
that the gauge-covariant derivative acts on the super- 
fluid order parameter, which is a complex scalar field 



the dynamics of the phase (p. It can also capture the 
strong-coupling phase if we incorporate the agents that 
destroy phase coherence: the vortices, windings of 2tt 
of the phase field. From the dual viewpoint, vortices are 
particles, that can condense just as well as bosons carP^^l, 
The vortex condensate corresponds to the state where the 
original variables ip have completely lost their meaning. 
In other words, the weak-coupling phase of the vortices 
is the strong-coupling phase of the original variables. 

It is useful to go over to dimensionless variables de- 
noted by a prime, 

Se = ^5*15, X = ax , T = T , yl„j = ——A . 

Cph e*a 

(10) 

We shall suppress the primes in the remainder. The 
dimensionless version of the action Eq. Q is 5e = 
/ drd^x C with, 

/: = -^(5f^-A,)^-^i^.V (11) 

Here the dimensionless coupling constants are, 
1 Ja 1 ha"-^ 

(12) 



Two quantities are of interest in the duality. The first 
is the current — -{dj^^^tp—A^^), related to the charged 



in the London limit (ps constant) is. 



/ps^e"^. Hence, the minimal coupling prescription supercurrent as 



^J^J^^ Then Eq. ilU can be 



(8) 

Here e* is the electric charge of one boson (one Cooper 
pair). To preserve gauge invariance, the temporal com- 
ponent of the gauge potential should have the same ve- 
locity factor as the derivative, and therefore we define 
AP*^ = (— i— K A). In addition the MaxweU action for 

M ^ Cph ' ' 

the dynamics of the electromagnetic field is included, 
which is governed by the speed of light c. Defining the 
electromagnetic field tensor F^^^ = {df^Ai, — di,Af^) where 
dfi = (l^r, V) and — (— iiy,V), the total action is, 



Se = / drd^x [~^Ja 



2-D 



e 



1 



1 

4/io ^''^ 
(9) 

This is the relativistic Ginzburg-Landau (Abelian-Higgs) 
model, where we have suppressed the potential terms 
~ a|\E'p-|-/3|\E'|'*, which are frozen out in the London 
limit of small amplitude fluctuations, as we will assume 
throughout. 



dualized by direct substitution into, 

^duai = Igwl - z«^(aPV - Af) - j-F'^,. (13) 

The second quantity of interest describes the Abrikosov 
vortices, which are singularities in the phase fleld p. For 
the remainder of this section we specialize to 2-fl di- 
mensions, for simplicity. Splitting the phase field into a 
smooth and a singular part, p = (/'smooth + <^sing, a vortex 
solution of winding number N satisfies, 

2nN = f dXf_, df,p> / dS^ tKv^id^d^,{psmoot\i + (Psing) 

JdS JS 

= j dS'„ eKutj.dudfj^psing, (14) 

by Stokes' theorem. The derivatives acting on a singular 
field do not commute. 



The superconductor is a Coulomb gas of 
vortices 



III. VORTEX DUALITY IN 2+lD 

In the previous section we derived the weak-coupling 
continuum limit of the Bose- Hubbard model in terms of 



On the smooth part we can perform integration by 
parts in Eq. (13), to obtain a term {d^^Wf^)ipsnwoth, and 
Vsmooth can be integrated out as a Lagrange multiplier for 
the contraint dj^^Wf^ — 0, the conservation of supercur- 
rent (continuity equation) in the superconductor. In 2-1-1 
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dimensions, this constraint can be explicitly enforced by 
expressing it as the curl of a dual gauge field, 



(15) 



This expression is invariant under the addition of the 
gradient of any smooth scalar field. 



(16) 



Substituting this definition in Eq. ( 13 ) leads to. 
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2+ld — ^yy^ti. 



+ e,,.df^b,Af~^^Fl,. (17) 

Here we have performed integration by parts on the sec- 
ond term, and we recognize the expression from Eq. ( 14 1. 
Therefore we define J^{x) = ei^^fj,d^dfj_(psing{x) as the 
vortex current. It is a one-form (vector) field because 
a vortex point particle traces out a worldline, with line 
element J^{x). Then from the coupling term b^J^ in 
Eq. ( jlT] ) we see that vortices interact by exchanging dual 
gauge particles b^. In other words: the 2-l-lD neutral su- 
perfluid (e* 0) is a Coulomb vacuum for the vortices 
with long-range interactions mediated via dual 'photons' 

In the charged superfluid, the current tu^ ~ ^fiVKd^^bK, 
also couples to the real electromagnetic photon A^j^, ren- 
dering the interaction between Abrikosov vortices short- 
ranged. However, for the formation of the vortex conden- 
sate as described below, it is unimportant whether or not 
it is coupled to electromagnetic field. This implies that 
the extension to the charged superfluid/superconductor 
is straightforward by choosing e* > 0. 



B. The Bose-Mott insulator is a dual 
superconductor 



condensate is a dual superconductor. We know what 
the phase transition implies for massless gauge fields: 
as they couple minimally to a condensate field $ = 
|$|e"^ they become massive due to the Anderson-Higgs 
mechanisrri^^'^^*^. We end up with the dual supercon- 
ductor in terms of the gauge field b^, 

+ -i6.)$p-^|$P-^|$|^ (18) 

where a and /3 are dual Ginzburg-Landau parameters. 
The gauge field 5„ obtains a Higgs mass and fur- 

thermore the longitudinal polarized photon now becomes 
a physical degree of freedom. This is the dual way of 
expressing that the single massless zero sound mode of 
the superfiuid turns into the two gapped modes (doublon 
and holon) in the Mott insulator^! Therefore the former 
is referred to as the Coulomb phase and the latter as the 
Higgs phase in terms of the vortex operators. 

In the charged superconductor, the original Anderson- 
Higgs-massive Goldstone mode (phase of the supercon- 
ducting condensate ^') dualizes into b^ and its gets con- 
tributions to its mass both from the electromagnetic field 
and from the vortex condensate. Furthermore the vortex 
condensate $ gives rise to an additional mode. In other 
words, the transverse polarization &t is the superconduc- 
tor sound mode, and the longitudinal polarization &l is 
the vortex condensate sound mode. They are gapped and 
degenerate. In addition there is still the electromagnetic 
field A^ which is also gapped. To avoid confusion by all 
the various gauge fields in this 'dual equation', it is of- 
ten useful to first regard the neutral limit e* — > and 
subsequently let the electromagnetic field enter weakly 
coupled to the current w^. 



The true power of the duality lies in the fact that the 
strong-coupling phase, i.e. wildly fiuctuating phase fields, 
can be described as an ordered state in terms of the vor- 
tices. Vortex-antivortex pairs can spontaneously emerge 
and annihilate in the form of closed spacetime loops. In 
the Coulomb phase (superfluid), such processes are heav- 
ily suppressed, as the coupling constant g acts as the line 
tension of such spacetime loops. 

As the coupling constant g decreases, the vortex world- 
lines grow in size and number, until at the critical point 
gc they span the whole system. At that point, vortices 
and antivortices can be created energetically for free. In 
other words, across the phase transition we find an or- 
dered state, the vortex condensate <&, out of which vor- 
tices can be pulled everywhere for free, just as Cooper 
pairs can be pulled out of the superconducting conden- 
sate. This is pictured in Fig. [2] 

Since the dual gauge fields b^ couple to the vortex 
condensate <I> just as the electromagnetic field cou- 
ples to the superconducting condensate the vortex 



IV. VORTEX DUALITY IN 3+lD 

The question is now how this generalizes to higher di- 
mensions. In the boson language, the superfluid/Mott 
insulator picture is unaltered. But vortices become ex- 
tended objects: Nielsen-Olesen (non-critical) strings in 
3+1 dimension^^H xhe Abrikosov vortex line traces out 
a worldsheet in spacetime, with surface element Jj^ = 
^KXvudud^ifsiug (Fig. [l|. The temporal components JjY 
are the density of the vortex line along I, while de- 
notes the motion in the direction k of the line along 
with continuity equations d^J^^ — for all A. 



A. The Coulomb phase 

The are two-form antisymmetric tensor current^^Sl^ 
and the dual gauge fields mediating the interaction be- 
tween vortices also become two-form fields bf^\. For 
the Coulomb phase (superfiuid), the generalization is 
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FIG. 2: Cartoon of the superconductor-Bose-Mott insulator transition in the dual vortex-language, for simplicity shown in 
2+1 dimensions. In the superconducting state of a system of bosons tunneling between potential wells, all the particles are 
completely delocalized (left). When the local repulsions become large enough at a density of one boson per well a "trafffic jam" 
sets in and the Bose-Mott insulator is formed (right). In the dual view one focuses on the physics of the vortices, the topological 
excitations of the superconducting state. In the superconductor these occur as virtual fuctuations of bound vortex-antivortex 
pairs, forming closed loops of worldlines in spacetime (red lines, middle left). However, at the quantum phase transition to the 
Mott insulator these loops "blow out"; the Mott insulator itself corresponds to a tangle of worldlines (middle right). Due to 
long-range interactions between vortices this just represents a relativistic superconductor (Higgs vacuum) of vortices, the Higgs 
mass being coincident with the Mott gap. Here we consider the situation in 3-1- ID where the vortices turn into strings. The 
Mott insulating state is now dually described as a 'foam of vortex strings' in spacetime, eventually behaving just as a three 
dimensional superconductor. 



freedonP3. It is the purely transversal component with 
K, and A each taking a transversal direction. For the elec- 
trically charged case this mode becomes gapped, just as 
those of the electromagnetic field do. Nevertheless, 
the superfluid dualizes into a gas of vortex worldsheets 
interacting via two-form gauge fields. Vortex-antivortex 
creation and annihilation events (quantum fluctuations 
in the superfluid) take the form of small closed world- 
sheet surfaces, suppressed by the large coupling constant 
9- 



B. The Higgs phase 

The dual condensate corresponds with a 'foam' formed 
by vortex strings fllling 3-l-lD spacetime. This 'stringy 
Higgs phase' is obviously somehow different from the 
conventional 'particle' Higgs phase of 2-l-lD. Surely the 
superfluid is ordered in terms of phase dynamics, and 
the Mott insulator corresponds to the completely phase- 
disordered state. A phase winding of 27r corresponds to 
the local formation of a vortex excitation. Therefore we 
expect that the Mott insulator is again a condensate of 
such vortex excitations. This cannot be dealt with us- 
ing standard field-theoretic techniques as in 2-f ID, where 
one can write down a quantum field theory of meandering 
vortex worldlines, the collective field of which takes the 
form of a Ginzburg-Landau scalar field^''. Now we should 
have a quantum field theory of vortex worldsheets: a 
string field theory. Such a theory is not yet available in 
closed form. 

However, at least for the mundane finite-energy vor- 
tices in condensed matter — as opposed to the coreless 



straightforward. The conservation of supercurrent 
djl^w^ = can be enforced by expressing it as the curl 
of this two-form gauge field, 

Wf, = et^unxdf'h^x. (19) 

Here e^uK,\ is the completely antisymmetric Levi-Civita 
tensor in four dimensions. This expression is invariant 
under the addition of the gradient of any smooth vector 
field, 

h^,x^h^x + d^,ex-dxe^. (20) 

Even though the transformation looks like the electro- 
magnetic field strength F^y = d^A^ — d^A^, it should 
not be confused with the actual field strength related to 



this dual gauge field, which is in Eq. ( 19 1. The gauge 
transformation itself has a redundancy, as another gauge 
transformation — sx + 9x1] would yield the exact same 
transformation, which is sometimes referred to as "gauge 
in the gauge" , and is important for the counting of de- 
grees of freedorrPS. 

The Lagrangian of the Coulomb phase Eq. (17) gen- 
eralizes to 

•^3+id — ^9{^tiVKxd^^biix)^ — b^xJ^x 

+ e.^^xd^b^xAf - ^^Fl,. (21) 

Again, the vortices have long-range interactions me- 
diated by the superfluid sound mode, now represented 
by 5k A- In the neutral case e* — 0, we retrieve the 
theory of a free and massless two-form gauge fleld in 
3-l-lD, which is known to have one propagating degree of 
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critical strings of string theory — the final result must be 
the Bose-Mott insulator. This insulator has two gapped 
doublon and holon modes regardless the dimensionality 
of the system. Hence, whatever the vortex string con- 
densate may be, it should add precisely one dynamic 
mode aside from the sound of the superfluid, and both 
modes should become gapped and degenerate. This is 
precisely the guiding principle we employed in our earlier 
work Ref. 1151 which we now briefly summarize. It turns 
out that earlier attemp ts at establishing a field theory 
of vortex stringsHEZED assumed that one can straight- 
forwardly generalize the minimal coupling construction 
for the phase of the Higgs field of second quantization 
{dfj_(p - b^) to the stringy case (i9^c,. - a^c^ - b^^). How- 
ever, this implies that one has to associate a vectorial 
phase field to the string condensate, which yields two 
longitudinal photons and three massive modes in total. 
Although this might be accurate for critical strings, it 
does not add up to the doublet of gapped modes of the 
Bose-Mott insulator. 

Therefore we reconsidered the status of the two-form 
gauge field bi^\. The single reason for introducing it in 



Eq. ( 19 ) was that the supercurrent is a conserved 
quantity. If we resubstitute this definition in the Higgs 
Lagrangian for 2+lD Eq. (18), we obtain. 



1 9 1, ,9 1 



1 

4^ 



fJ.V 1 



(22) 



The second term explains the Anderson-Higgs mecha- 
nism in the sense that supercurrent can no longer be cre- 
ated for free (the modes are gapped/massive), but it does 
not explicitly demonstrate where the additional degree of 
freedom, the 'longitudinal photon', originates. 

We need to realize that a vortex is a source or sink of 
supercurrent. Therefore, in the vortex condensate where 
vortices can be created for free at every point in space, 
the conservation of supercurrent is violated. A more pre- 
cise statement is that there is a superposition of having 0, 
1, or any number of vortices at any point, such that cor- 
relations of the phase field vanish completely, a notion 
we explored further in Ref. [221 Hence, the constraint 
dj^w^ = is removed, which liberates the longitudinal 
component of the current as a physical degree of freedom. 

Therefore the Lagrangian Eq. (22) is valid in any 
dimension. The vortex condensate amounts to the ap- 
pearence of the second term ^ the Higgs mass / 

condensate density / Mott gap. Concurrently, the super- 
current is no longer conserved, and the longitudinal com- 
ponent of the supercurrent enters as a physical degree of 
freedom, leading to two gapped modes in the Bose-Mott 
insulator, in any dimension. 

The electromagnetic field couples as always to the elec- 
tric current — e*w^. 

Summarizing, the 3-l-lD Bose-Mott insulator is again a 
dual superconductor, albeit of a special kind where two- 
form gauge fields take the role of Higgsed photons. Nev- 
ertheless, the dual order parameter <I> instigates a dual 
Meissner effect by causing electric current to decay ex- 



ponentially resulting in the insulating behavior. It also 
immediately suggests that the vortex condensate has vor- 
tices J^^ = eK\fj,udf^df;^(j) of its own, which are lines of 
quantized electric current just as Abrikosov vortices are 
lines of quantized magnetic flux. This we will investigate 
in further detail below. 



C. Dual vortices and the dual gauge field 



Still, the question remains how Eq. ([22j) can be ex- 
pressed in terms of the two-form gauge field b^^x- This 
issue is particularly important considering vortices in the 
dual condensate. What are the singularities in the phase 
field (j) of the stringy vortex condensate order parameter 
$ = l^le'"^? 

All along the problem is how to match the gradient of 
this phase field to the two-form gauge field. What is the 
form of the minimal coupling analogous to {d^ — ib^)^ of 
Eq. pSl): 



{d^ ~ i ??? 6,a)$ ? 



(23) 



Earlier woriy^liZiasi proposed a vectorial phase field 
(j) c^, but this implied a too large number of degrees of 
freedom as we already mentioned above. In our preceding 
workl^ we proposed. 



id. 



ie^lkA^KA)*, 



(24) 



where the three free indices in e^n^A take values orthog- 
onal to the four-momentum 15^ — >■ only. Essen- 
tially this corresponds to the generalized Lorenz gauge 
fix d^bnx = VA. This does lead to the correct form Eq. 
( p2| ). However, the relation between the phase field 
and the dual gauge field 6kA is obscured, because cross- 
terms are explicitly eliminated by this gauge fix. The 
same thing happens in the Ginzburg-Landau equations 
in the Lorenz gauge fix, the phase and the photon fields 
become decoupled. Any inquiry into the vortex excita- 
tions, which are singularities in the phase field, cannot 
rely upon such gauge-fixed expressions. 

It turns out that the issue of worldsheets of current 
lines is actually a surprisingly intricate affair. One would 
think that for instance a simple current-carrying cop- 
per wire cannot hold any macroscopic secrets, but the 
equations of motion when such a wire is regarded as a 
worldsheet in spacetime have never been established. In 
an earlier work'^*^' we performed such an investigation for 
the case of Abrikosov vortices in superconductors. Here 
everything fits together neatly. Start out with the rela- 
tivistic equation of motion derived from Eq. (|9|, 



(25) 



where we suppressed all dimensionful constants. Now act 
on this expression with e^xp^dp to find. 



d ^nXfii^-Fpiy ~i~ ^nXpi'-Ffiiy — '^kX 



(26) 
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Here is the density of an Abrikosov vortex line in 
direction I, and is the motion in direction k oi a 
such a vortex along I. By neglecting the first term 
which denotes the decay of the electromagnetic fields 
away from the vortex due to the Meissner effect, we see 
that a vortex line J^Y induces or couples to magnetic field 
Bi = £kmnFmn, and motiou of the vortex induces elec- 
tric field e-kitnFtn = eunEn, corresponding to the 
relation v x B = E, well known in vortex dynamics. 

The reason why this works so well is that the two- 
form vortex source Jj^ couples directly to the two-form 
electromagnetic tensor i^^^. Conversely, the electric cur- 
rent is a vector quantity. So we face a problem similar 
to Eq. ( 23 ) , matching quantities of differrent geometrical 



nature. A closed-form relativistic expression is not found, 
but the physics at play can nevertheless be accurately 
comprehended. This is illustrated in Fig. [3] imagine a 
volume element of current at position x. This current el- 
ement cannot distinguish between moving along a static 
line (wire) in a certain direction I, or being dragged along 
a vortex line along k but moving in direction I. There- 
fore, the current in direction I gets contributions both 
from static lines pointing along I, denoted by JT^Y, and 
from lines pointing in different directions k moving in di- 
rection /, denoted by J^. Our task is to, given a vortex 
worldsheet J^^, derive the resulting electric current w^. 

It is instructive to first observe the 3-l-OD limit of static 
vortex lines carrying stationary current. This is in fact 
also the starting point for discussions of Abrikosov vor- 
tices. In the static limit the only vortex components are 
= = etimnd^d^(l), and the only dual gauge field 



components are bi = hti- The Lagrangian is Eq. (18) 
where however all indices denote spatial dimensions. It 
contains the standard form of the minimal coupling, and 
hence the equation of motion is readily obtained by vari- 
ation with respect to hi, 



(27) 



Here Bi = ei„md„Am is the magnetic field. Acting on 
this expression with ersids gives, 

g{drd,^-5r^d'')Wra + \^?Wr+er,ldsBi = . (28) 

Here we clearly see that a static vortex line sources 
electric current Wr in the parallel direction. This current 
decays exponentially in the insulator. The current also 
sources the curl of magnetic field, as in the Ampere- 
Maxwell equation but we will see later that this effect is 
very weak. 

This is the situation of Fig. 



3(a) It is clear that by in- 



corporating the time dimension, the vortex line traces out 
a worldsheet, the four-current has the charge density as 
temporal component, and hence Fig. 3(c) describes pre- 
cisely the same situtation as Fig. 3(a) The worldsheet 



surface element is spanned by the temporal component 
Wt — cp and the spatial current w;, such that the four- 
current is 'the diagonal' of the worldsheet. Now consider 
the situation of Fig. |3(d)[ Here the worldsheet is 'tilted' 



and therefore contains, next to the density J^^, a com- 
ponent This surface element is spanned by the two 
current components and w^- 

These notions are most conveniently expressed by us- 
ing differential forms (see any text book on differential 
geometry, e.g. Ref. 26). Here dx^ represents a line el- 
ement in direction x^. A vector field a(a:) is expanded 
in components as a(a;) = fj, o. ii{x)(ix'^ = a^{x)Ax^, and 
integration follows directly Ja(a;) = J a^j^{x)dx^^. A vec- 
tor is in this context called a 1-form. We can combine 
multiple 1-forms into n-form by using the wedge prod- 
uct A, which is the antisymmetrized tensor product ®. 
Thus da:^ A dx" = da;^ ® dx" - dx" (g) da;''. These defini- 
tions show that a 1-form describes a line element and a 
2-form describes a surface element, and this language is 
therefore perfectly suited to describe our worldsheets. 

Strictly speaking, the vortex worldsheet is a d — 2-form 
and the current is a d — 1 form (see Ref. [30|) . But we im- 
plicitly use the so-called Hodge duals of these quantities, 
turning them into 2- and 1-forms respectively, defined by 
= da:" A dx^ and w = tu^da''. Based on the 
considerations laid out above, we propose that 

the vortex worldsheet and the current = [cp, 



are related as 



J\ da'' A dx^ 



1 

cp 



w^dx" A uiAda 



(29) 



The factor — is necessary to obtain the correct dimen- 
sions, such that — -^wtWi — wi is a charge per time. 



agreeing with Eq. (28 1. In appendix [Bl it is shown that 



this relation obeys the desired behavior under Lorentz 
transformations. We must conclude that, even though 
there is no closed form for the vortex equation of motion 
like there is for Abrikosov vortices Eq. (26), the physical 
content is completely clear. Once the Mott condensate 
has formed, near the quantum phase transition it allows 
for vortex excitations which are lines of quantized current 
that obey Eq. (pgl). 



Let us emphasize again that we have assumed here the 
strong type-II limit, corresponding to the London limit 
where fluctuations in the condensate amplitudes jvl/l, |$| 
are suppressed. Phenomenologically this means that the 
penetration depth A is much larger than the coherence 
length ^. It is well understood how to interpret the dif- 
ference between type-I and type-II behavior in the duality 
context. In 2+lD (or 3D classical)!^ the type-I regime is 
associated with effective attractive interactions between 
the vortex particles, triggering a Van der Waals-type first 
order liquid-gas transition, translating into the first order 
transitions of the type-I state. The same logic of course 
applies to the 3+lD context where net attractions be- 
tween vortex worldsheets would have a similar effect. We 
just assume that the interactions are repulsive, prohibit- 
ing the clumping of vortex matter, such that the systems 
exhibit the continuous phase transition of the 3+1D/4D 
XY-model, equivalent to the type-II regime. 
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(a)static vortex line 



(b)dynamic vortex 
pancake 



(c)static vortex 
worldsheet 



{d)dynamic vortex 
worlsheet 



FIG. 3: (a) Static vortex line in the xy-plane; the current flows through the line, (b) Vortex pancake moving in time (blue). 
The associated current in the spatial direction is shown in red. |(c)| Static vortex line in the a;2-plane moving straight up in time, 
(d) A vortex line in the z-direction moving in the a;-direction through time. The last two worldsheet configurations correspond 



to the same electromagnetic current (red). 



PHENOMENOLOGY OF THE BOSE-MOTT 
INSULATOR 



A. Maxwell equations 



In this section we derive observable quantities of the 
Bose-Mott insulator and its vortices. The lack of a com- 
plete expression for the dual gauge fields is not as large 
an obstacle as one might think. The remainder of the pa- 
per will focus on the 3-l-OD limit of static vortex lines of 
stationary current, which will prove interesting enough. 
The complete dimensionful Lagrangian corresponding to 
Eqs. (pl and (pSl) is, 



The last equation Eq. (33) is clearly the inhomoge- 



neous Maxwell equations for a source term — j-Wm- 
These equations carry over from the superconductor, and 
do not pertain as such to the Mott insulating state. The 
insulating behavior is due to the expulsion of the elec- 
tric current, which is represented by the term ~ I'&p. 
Therefore, Eq. (331 is just the vacuum contribution to 



dynamic electric and magnetic fields generated by a cur- 
rent source. 



he, 



2a^ ^ ncnh 



-,D-2 



2J 



e* 1 
n 4^0 



ph 

|$|2_ ^|$|4 



(30) 



Here |<i>p is the dimensionless condensate density. In 
what follows we will consider D = 3 only. 



We derive the equations of motions by varying Eq. ( 30 ) 
with respect to bk and The equations of motion 
are. 



ch 

$-5$-/3|$|2$ = 0, (31) 
a |$p a 1 

' ^^knm^n'^m ~l~ \^k^ 11 ^k) ~" ~^^k7nn-^ mm 

J a /iCph ^ 

(32) 

1 e* 
Mo n 



Here we have substituted the definition Wm — ^mnkdnbk- 
We are now ready to discuss the physical content of these 
equations. Note that the last two equations reduce to 
the equations of motion associated with the standard 
Ginzburg-Landau superconductor in the limit |$| — >■ 0. 



B. Penetration depth 

The dual penetration depth Am sets the length scale 
up to which an electric current penetrates in the Mott 
insulating region. To identify it we act on Eq. (32) with 
trskds- Contracting repeated indices, and using drWr = 
0, we find in the London limit of the dual superconductor 
l$l = 



«a2 * 



he. 



-Wr 



'ph 



e 

—d F 

h 



$2 



a 



(34) 



Here we recognize the Mott vortex current 
^rskdsdk4'- The interpretation of this equation is as fol- 
lows: a supercurrent Wr can be generated by a vortex 
source . This current is 'dual Meissner screened' by 
the Mott condensate $00 as witnessed by the second 
term, but there is also some electromagnetic screening 
from the 'backreaction' of the induced electromagnetic 
field. In order to see this, we substitute Eq. (33) in this 



equation. In the absence of vortex sources, this leads to, 



pa tEM 



EM 



$2 



tEM 



he-pY, 

hps 

Cphm 



h-^ 



-J, 



EM 



1 



tEM ^ T 



EM 



= 0. 



(35) 



Here we substituted a/ J = m*/h Ps [see Eqs. (12 1 and 
([6])], and used the definition of the London penetration 
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depth Al = jjiQe*"^ ps/m* . We find indeed two contribu- 
tions to expulsion of electric current. The first ^ <i>^ 
is due to the Mott insulator, and the second remembers 
that the system originated from a superconductor. This 
is actually rather odd: the Meissner screening is due to 
the fact that the superconductor wants to expel the mag- 
netic field, which is not true for the Mott insulator. Again 
one must consider that the insulator is a phase-disordered 
superconductor, and that on (very) short length scales 
the local boson superconductor is retrieved. Let us make 
a crude estimate of the relative strengths of the screening, 
by inserting the numerical values. 



/^o = 47r.l0"^ w lO^^N/A^ 



IQ-i^C, 



h « lO^-'^Js, Cph 



1 

300' 



lOV/s 



we find that the relative strengths are. 



Mott 



4^2 



Meissner /ige* Cph/^ 



10*$ 



(36) 



(37) 



Here $^ is dimensionless, but as the order parameter 
of the Mott condensate it should be large. Therefore 
the expulsion of electric current due to the Mott term is 
several orders of magnitude stronger than the Meissner 
screening, and for all purposes the latter may be ignored, 
also eliminating our interpretation problem. 

Hence the "Mott proximity depth" for electric current 

is Am = \I * Ps^L- It depends on a number of mate- 



Cph 

rial parameters. We encounter the ubiquitous combina- 
tion Ps^lo: which is the product of the superconducting 
order parameter and the Mott order parameter. At first, 
one may think that they should be mutually exclusive, as 
one has either superconducting order or Mott insulating 
order. However one must realize that the Mott insula- 
tor is formed from repelling Cooper pairs: the larger the 
number of Cooper pairs, as quantified by the superfluid 
density ps, the stronger the electromagnetic effects such 
as screening will be. It is just that signals the exis- 
tence of the Mott state, whereas the combination Ps^^ 
is the appropriate dual Higgs mass of this insulator. 



C. Coherence length 



Upon rescaling the dual order parameter <i> in Eq. (31 ) 

by dividing it by its equilibrium value $oo = \J~^^ such 

that $ = $oo'l'', and taking 6^ = as is the case deep 
within the Mott insulator, this equation reduces to. 



Landau parameter \a\. Fluctuations in the order parame- 
ter will take place within a typical length ^. In the strong 
type-II limit Am > 



D. Current quantization 

We are now in the position to present our central result: 
the quantization of vortex lines of electric supercurrent. 
Eq. ( 32 1 is similar to the regular Ginzburg-Landau equa- 



tion. Consider a closed contour over which the change of 
the phase is a multiple of 27r, 



dx'' dk(, 



2tiN. 



(39) 



dS 



We are free to choose this contour deep within the Mott 
insulator far away from the vortex line, such that the 
electric current is suppressed Wm = 0. Assume that there 
is no external electromagnetic field F^^^ = 0, while the 



induced field is very small as argued in Eq. (37). Under 
these conditions Eq. (32) reduces to. 



he- 



ph 



(40) 



Taking the line integral of this equation as in Eq. ( 39 ) , in 



voking Stokes' theorem on the right-hand side, one finds. 



he. 



ph 



2itN 



he 



ph 



dS 



dx" bk 



ds 



— / dSrn emnldnbl — / d5„ 

Js Js 



(41) 



The right-hand side is the fiux of current w,„ through the 
surface S. Since the current is expelled from the Mott 
insulator, this current flows through the vortex line. For 
the electric current I which is the fiux of the current den- 



sity J„ 



-Wm, this implies the quantization condition. 



In - 



.^2^7V= — Vc/J(27r) 
a $0 



■N. 



(42) 



Here $o — h/e* is the (magnetic) flux quantum and 
we have substituted the microscopic parameters VUJ — 
hcph/ci from Sec. |ll] Following similar arguments as for 
the Abrikosov lattice, higher winding number {N > 1) 
vortices are energetically unfavorable, and split up into 
multiple iV = 1 lines. Each such vortex line carries an 
electric current, for typical values Cph ~ 10® m/s and 



10- 



-10 



m, of. 



Cph 



27r « 10-^A. 



(43) 



^(a„)2$' + $', 
\a\ 



0. 



(38) 



This shows that we can define a dual coherence length 



as ^ 



depending on the details of the symme- 



try breaking through the precise value of the Ginzburg- 



Different from the magnetic flux quantum, this cur- 
rent quantum is not exclusively expressed in fundamental 
constants for the obvious reason that the unit of current 
(Ampere) cannot be composed this way. However, we 
do observe that the current quantum is inversely propor- 
tional to the flux quantum in Eq. ( 42 1 as expected by 
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duality. However, one needs a quantity with dimension of 
energy (y/UJ) to convert the reciprocal of the flux quan- 
tum into a quantum carrying the dimension of charge per 
time. This then conspires into the combination of su- 
perconducting phase velocity and lattice constant of Eq. 
(42). Although depending on material specifics, Cph/a 



is expected to be a fixed quantity in the proximity of 
the superconductor-Mott insulator quantum phase tran- 
sition. 



VI. PHASE DIAGRAM 

The results in the above lead us to propose a general 
phase diagram for the charged Bose-Hubbard model in 
3-l-lD (Fig. |4|. The control parameters are i) the (quan- 
tum) coupling constant g ~ \/U/ J ; ii) temperature; and 
iii) an applied magnetic field for the superconductor, or 
an applied electric current for the Bose-Mott insulator. 

The superconducting side is well known. At [/ ^ 
we have the familiar _ff-T-diagram for type-II supercon- 
ductors (we always assume local pairs, ^ <^ Al). At low 
applied field the Meissner phase is formed with complete 
field expulsion. Above Hd ^ <&o an Abrikosov vortex 
lattice forms, until the superconductivity is destroyed at 

Upon approaching the superconductor-Mott insula- 
tor transition by increasing the coupling constant g, the 
renormalized phase stiffness is decreasing and the super- 
conducting Tc vanishes right at the quantum phase tran- 
sition. 

Contrary to standard expectation that on the Mott 
insulating side of this quantum phase transition a fea- 
tureless state is found, the first surprising ramification of 
our duality is that in 3+lD there is actually a thermal 
phase transition in the Mott state with a finite critical 
temperature that is rising (at zero applied current) upon 
moving away from the QPT. As usual, the thermal phase 
transition occurs on both sides of the QPT well inside the 
renormalized classical regime; upon raising temperature 
one will re-enter the quantum critical 'wedge'. 

Duality dictates that the magnetic field in the super- 
conductor turns into electric current for the Bose-Mott 
insulator. Hence, the control parameter equivalent to 
the magnetic field of the superconductor is applied cur- 
rent for the Bose-Mott insulating side. Near the quantum 
phase transition where duality is expected to be univer- 
sally valid, this suggests a dual Meissner phase (insula- 
tor) below the current quantum /q, while at higher cur- 
rent a dual Abrikosov lattice of quantized current lines 
will occur. When current gets too high, above the up- 
per critical current Ic2, the Bose-Mott insulator will be 
completely destroyed. The resulting state is again a su- 
perconductor, since the dual of the Bose-Mott insulator 
is a superconductor and not a metal. Of course, if the ap- 
plied current exceeds the pair-breaking current density of 
the superconductor, superconductivity is destroyed, and 
the Bose-Hubbard model no longer makes sense. 



Upon increasing U /J further, eventually one ends up 
in the simple strong-coupling limit of the Bose-Mott in- 
sulator, with the nearly completely localized bosons sub- 
jected to short range virtual fluctuations. This is indeed 
a featureless state, one that sets the common intuition. 
How to accommodate the rich physics we claim near the 
quantum phase transition? This is subtle, but in fact 
quite simple. In the language of the dual superconduc- 
tor, the dual penetration depth shrinks to the lattice con- 
stant: the effect is that our current-carrying vortices 'fall 
through the lattice constant'. Their core energy exceeds 
the UV cut-off and the the "type-I Bose-Mott insulator" 
can be formed, with a Tc that has 'disappeared above the 
cut-off'. In this limit the state turns featureless again, 
just governed by the thermal excitations of the 'massive 
photons' (doublons/holons) of the dual Higgs condensate. 

The strong coupling is therefore part of the reason for 
the misleading intuition that the Bose-Mott insulator has 
to be thermodynamically featureless. Another flaw in 
this regard is associated with approaching it from the 
other side: it is well understood that the quantum phase 
transition in 3-f ID is at the upper critical dimension, as 
it is effectively XY in 4D. The critical regime is therefore 
governed by mean-field and dominated by amplitude fluc- 
tuations. How then can these vortices play such a central 
role? 

This is just a confusion based on overestimation of 
the influence of universality class away from the crit- 
ical regime. A central result of the renormalization 
group/critical theory is that the "soft spin" (1^9^) and 
"hard-spin" (sigma model, our Bose-Hubbard model) 
share the same, universal critical regime which is surely 
of the mean-field kind. However, duality is actually re- 
lating the (stable) fixed points on both sides of the phase 
transition: our dual superconductor becomes discernible 
only well below the cross-over to the quantum critical 
regime. The quantum critical regime is itself governed 
by the mean-field dynamics of 4D XY. 

Surely, starting with microscopic circumstances that 
directly coarse grain in a i^s^-theory, there will be no in- 
teresting physics on the insulating side: this physics plays 
no role in a real BCS superconductor (such as aluminum) 
since amplitude fluctuations dominate on all scales. How- 
ever, starting out with strongly bound, hard-core bosons, 
duality cannot be avoided. All that matters is that the 
dual penetration depth be large compared to the lat- 
tice constant near the quantum phase transition. This 
length scale is coincident with the typical distance be- 
tween free vortices and as long as this is large compared 
to the lattice constant the Mott insulator must be the 
dual vortex superconductor. With regard to the critical 
regime one just learns that the proliferation of the vortex 
strings is described by a mean-field regime, regardless of 
whether these strings carry magnetic fiuxes (supercon- 
ducting side) or electric currents (Mott side): it remains 
an interesting exercise to find out why the se cr itical the- 
ories match in one common critical regimd^'^. 

Summarizing, the type-II Bose-Mott insulator is a dual 
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FIG. 4: Topology of the phase diagram of type-II Bose-Mott insulators in three space dimensions. The horizontal axis 
represents strength of quantum phase fluctuations that tune from superconducting (SC) to Bose-Mott insulating (BMI) order 
via a quantum phase transition. As function of temperature (vertical axis) the "wedge" of quantum critical fluid (QC) emerges 
from the quantum critical point, while on either side thermal phase transitions occur: obviously to a superconducting state, 
but duality also predicts a thermal transition on the insulating side. An external magnetic field is first expelled by the 
superconductor (Meissner), while above Hd it penetrates as an Abrikosov lattice of quantized flux lines. Near the quantum 
phase transition a similar response is observed in the Mott insulator, now as function of applied electric current. First the 
current is expelled, just an insulator. However, above a 'lower critical current' Id current will penetrate as an 'Abrikosov lattice' 
of lines carrying a quantized current flux: the type-II Mott insulator. Upon increasing the local repulsions one eventually enters 
the regime of nearly completely localized bosons. Here the type-II Mott insulator disappears because the relevant length scales 
have shrunk beyond the lattice constant: the reason the type-II phase has hitherto been overlooked. 



type-II superconductor. The dual of the magnetic field 
of the superconductor, is electric current for the Bose- 
Mott insulator. To probe a superconductor, one applies a 
magnetic field from outside. Outside the superconductor 
there is a medium which supports a magnetic field with 
magnetic permeability fi. It does not matter whether this 
is a vacuum with permeability fiQ or a dielectric with an- 
other value. Even a metal would do when considering a 
static magnetic field. As such, there are only two inequiv- 
alent phases with regard to magnetic field: Meissner and 
not-Meissner. 

For electric current, the situation is different. The 
superconductor is to the Bose-Mott insulator what the 
(Maxwell) vacuum is to the superconductor. The su- 
perconductor carries the current applied from outside to 
the Bose-Mott insulator. This begs following the ques- 
tion: What is the Maxwell vacuum to the Bose-Mott 
insulator? The vacuum has no charge carriers, and does 
not support a current. Also a metal is different, since 
a dissipative current of quasiparticles (fermions) may be 
different from a supercurrent made out of Cooper pairs. 
Therefore we anticipate that the situation for the type-II 
Bose-Mott insulator is richer than that for the supercon- 
ductor. This was already mentioned in Table |l] and is 
also exploited in the experimental setups of the next sec- 
tion. 



VII. PROPOSED EXPERIMENTS 

It is not an accident that this type-II Bose-Mott insu- 
lator was never seen in the laboratory. The best model 
systems are either hard to realize in three dimensions 
(Josephson networks) or it is unclear how to impose 
external currents (cold atoms), while in conventional 
condensed matter systems it is uncertain whether such 
physics is at work at all. In general one should focus on 
systems with large phase fiuctuations. Obvious candi- 
dates are strongly underdoped cuprate superconductors. 
Here the elusive pseudogap phase is by many conjectured 
to consist of so-called preformed Cooper pairs, which are 
bosons. Thus at some high temperature T* the pairs bind 
into bosons, and only at some lower temperature Tc phase 
coherence sets in, leading to superconductivity^'^^ 
Hence, the quantum phase transition from superconduc- 
tor to Bose-Mott insulator would be precisely of the 
vortex-proliferation kind discussed here. Experimental 
support comes from the Nernst effeclP^, diamagnetic 
behavioi'^ and spectroscopy^^'. One can also wonder 
whether the "giant proximity effect" associate d wit h a 
lOOnm thick underdoped cuprate barrier laye-^^^ has 
dealings with type-II Bose-Mott behavior (see Appendix 
[a] ) . If the pseudogap indeed consists of phase-incoherent 
local bosons, the type-II Bose-Mott insulator should be 
found close to the quantum phase transition and at low 
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temperatures. Conversely, if the current line lattice is 
found in the underdoped cuprates, it would constitute 
indirect but convincing evidence of the existence of pre- 
formed Cooper pairs. Another candidate may be the 
so-called disordered or amorphous superconductors, see 
for instance Ref. |39l 

Even though the imagination of experimentalists 
would be more qualified to devise the most suitable setup, 
we wish here to sketch some ideas to probe the current 
line lattice. In Fig. [5] we illustrate a number of possi- 
ble devices to measure the type-II Bose-Mott insulator, 
all revolving around the basic ingredient of imposing an 
external current on the Mott insulator. 

Setups (a)-(c) involve the classic junction type. One 
tries to force current through a slab of type-II Bose-Mott 
insulator. The first experiment (a) is analogous to the 
type-II superconductor. There, magnetic fields lines pen- 
etrate from the outside to form the first Abrikosov vortex. 
The current at first flows around the type-II Bose-Mott 
insulator, and when the current is large enough (higher 
than the current quantum Iq) it will 'trickle in' from the 
outside to form the first current line. The thickness of the 
superconductor surrounding the type-II Bose-Mott insu- 
lator should be smaller than the penetration depth Al, 
since in a superconductor current flows near the edges. 
Experiment (b) is similar, but there are no supercon- 
ducting walls. This presumes that current lines will form 
internally, not coming in from the outside. Below the 
lower critical current Id no current will flow at all, mak- 
ing for a stronger signal. Experiment (c) has normal leads 
instead of superconductors; nevertheless possibly any ap- 
plied current suffices to create vortex lines of supercur- 
rent. In these experiments, the slab of type-II Bose-Mott 
insulator should be thick, to preclude any normal Joseph- 
son effect. In all these cases, the measured signal would 
be a sudden jump in the current when a current vortex 
line is formed. 

Experiment (d) connects with so-called Joscphson vor- 
tices in superconductors. In a Joscphson junction, a nar- 
row barrier in between two superconductors, vortices of 
quantized magnetic field along the junction can form un- 
der applied field, which behave like Abrikosov vortices 
except that they do not have a normal core. Here a 
small barrier of superconductor is sandwiched in between 
two layers of type-II Bose-Mott insulator. The current 
through the superconductor would be quantized if the 
dual phase coherence of the type-II Bose-Mott insulator 
imposes on the narrow barrier. 

At first sight the SQUID setup (e) looks particularly 
promising. One would like to impose a current bias, and 
in the junction-type experiments above, one actually ap- 
plies a potential bias. One of the consequences could be 
that, since in general (Bose-)Mott insulators are rather 
poor insulators, leak currents may spoil the signal. In 
setup (e) a ring of superconductor is interrupted by a 
thick layer of type-II Bose-Mott insulator; this is a typical 
superconducting quantum interference device (SQUID), 
except that the barrier is intentionally very thick to pre- 



clude the normal Joscphson effect. Applying a magnetic 
field through the SQUID loop will cause a phase differ- 
ence across the insulating barrier. This phase difference 
will not induce a normal dissipative current, but if this 
difference is large enough, a vortex line of supcrcurrent 
may form. Then a current will flow through the loop 
with magnitude of one current quantum. Measuring the 
magnetic field through the loop (for instance with a sec- 
ond SQUID), one would see a sudden drop when this 
current starts to fiow. Increasing applied magnetic fur- 
ther would induce more and more current lines. The 
current quantum in a high-Tc material is estimated to 
be quite large [Iq ~ 10~^A, see Eq. (43)]. However, we 
would not be surprised when the pinning of these current 
lines would turn out to be very strong, given for instance 
the strong spatial inhomogeneity of the superconducting 
order in this regime as observed by scanning tunneling 
spectroscopy®^. This might cause substantial "current 
flux penetration" difficulties, in analogy with the com- 
plications that are we ll doc umented in the context of the 
usual vortex dynamic^^^'^. As for the Abrikosov lattice, 
the most direct way to probe the type-II Bose-Mott in- 
sulator would be the analog of decoration experiments, 
directly imaging the current lines. Scanning tunneling 
spectroscopy is here an option with the caveat that the 
size of the current line is set by the Mott proximity depth 
which can be quite large near the QPT. Alternatives are 
microwave impedance or low energy electron microscopy 
measurements. 



VIII. CONCLUSIONS 

We have demonstrated that Bose-Mott insulators 
which are close to the quantum phase transition to the su- 
perconductor exhibit a much richer physics than the intu- 
ition that follows from the strong-coupling limit would in- 
dicate. According to a precise quantum field-theoretical 
duality, its physics should be a mirror image of the rich 
physics of normal superconductors. The highlight is our 
prediction of the existence of an analog of the type-II 
phase, where now an Abrikosov lattice is formed of topo- 
logical vortex lines that carry a quantized flux of super- 
current. Since these Mott insulators have to be three 
dimensional, while the type-II phase is induced by cur- 
rents that are imposed from the outside, it is not straight- 
forward to see these effects in model systems that are 
designed to represent the Bose-Hubbard problem in the 
laboratory (cold atoms, Joscphson networks). However, 
the type-II effect can be exploited to find out whether 
such physics is indeed at work in underdoped high-Tc 
superconductors. It is imaginable that there is a range 
in dopings near the quantum phase transition where the 
low-Tc and low superfiuid density superconductors will 
turn out to be Mott insulators camouflaged as supercon- 
ductors due to a glassy current line network induced by 
the measurement fields. 
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FIG. 5: Various devices to measure the type-II Bose-Mott insulator. Indicated are the type-II Bose-Mott insulator (MI) in 
yellow, the superconductor in red and the Maxwell vacuum in blue, with the current lines as dashed black lines. The circle 
and arrow represent a current source, (a) The precise dual of the type-II superconductor: the outside superconductor carrying 
current acts like the Maxwell vacuum carrying a magnetic field, which is then penetrating the type-II MI in the form of the 
quantized flux lines. A complication is that the currents in the outside superconductor decay over the London penetration 
depth, (b) A Josephson junction using the MI as barrier. For a very small bias the barrier will be insulating, to become 
completely transparent above Id- (c) The medium imposing the current need not be a superconductor: a metal works as well. 
A simple "capacitor" will short circuit at Id. (d) Equivalent of Josephson vortices where the quantized vortex line does not 
form inside MI but within a narrow junction layer of SC. (e) SQUID setup in which current bias is increased in very small steps 
by a perpendicular magnetic field (circle with dot). Current will not flow until the first vortex is formed. 
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Appendix A: Giant proximity effect 

In the regular Josephson effect, a supercurrent can flow 
between two superconductors, even if there is a spatial 
gap or barrier in between them. Because the supercon- 
ducting order parameters extend outside the supercon- 
ductor, if the barrier is narrow enough that the two or- 
der parameters overlap, the supercurrent is supported 
even within the barrier. The order parameters fall off 
exponentially with typical tunneling length scale that is 
microscopic. If the region in between is a good metal, 
instead the normal metal coherence length is the ap- 
propriate length scale. 

Now, in type-11 superconductors and especially in the 
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FIG. 6: The geometry of the experiment of Ref. [37|on the gi- 
ant proximity effect. Between two layers of Lai.85Sro.i5Cu04 
(Tc = 40K) is one layer of La2Cu04+i (Tc = 20K) that can be 
up to more than 100 A thick. The yellow spots indicate gold 
leads over which a voltage bias can be installed. At T = 30K 
supercurrent was observed to flow (indicated by the dashed 
profiles) even through these thick barriers that greatly exceed 
the coherence length of only several A. In Ref. |38]a similar ge- 
ometry was used with leads of Lai.84Sro.i5Cu04 (Tc « 32K) 
and the layer in between of Lai.94Sro.o6Cu04 (Tc < 5K) that 
is 46nm thick. Muon spin resonance measurements at decreas- 
ing temperatures showed an increasing diamagnetic response 
to an applied field of 9.5 mT. 



high-Tc cuprates, the tunneling length is very short (of 
the order of several A). It came therefore as a sur- 
prise that placing an underdoped, non-superconducting 
cuprate layer between two cuprate superconducting leads 
(see Fig. [6]), a current was observed to flow even if the 
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barrier was very wide, up to hundreds of Af. A review of 
this "giant proximity effect" can be found in Ref. 143l It 
was speculated that the superconductor would Joseph- 
son couple between impurities throughout the barriers, 
but the very clean samples of Ref. .37, seemed to be the 
definitive answer that really a new phenomenon comes 
into play. Several explanations were proposed^^"^, all 
of which have in common that superconducting order is 
induced homogeneously within the barrier region. 

When the barrier region is a type-II Bose-Mott insu- 
lator as proposed in this work, the giant proximity effect 
will be established in a qualitatively very different man- 
ner. The supercurrent would penetrate as vortex lines 
of electric current, leading to very inhomogeneous con- 
ductivity as the vortices form a current line lattice. This 
explains very simply how the Josephson effect is "giant" : 
the energy cost of the formation of a vortex line is linear 
in its length, while the regular Josephson effect is lim- 
ited by the exponential fall-off of the order parameter. 
A clear test to confirm this prediction would be to ob- 
serve the spatial variation of the conductivity once the 
supercurrent is flowing. 

On a related note, very recently it was observed in 
sandwich samples very similar to those considered in 
the previous paragraphs, that the barrier region through 
which (super)current flows undergoes a Meissner effect: 
applied magnetic field perpendicular to the c-axis of the 
jmiction, is expellecj^by a diamagnetic response. There- 
fore we may regard this region to be superconducting. 
This agrees nicely with our proposed phase diagram, 
where an applied current above the higher critical cur- 
rent /c2 destroys the Mott insulator and at low enough 
temperatures drives the the system back to the super- 
conducting state. The applied field would induce a coun- 
tercurrent in the superconducting leads that causes the 
type-II Bose-Mott insulator in between to permit current 
lines, or is even pushed above its upper critical current Irji 
to become completely superconducting. The authors of 
this work claim that electric current must also be flowing 
in the a&-plane in which the cuprate layer lies, in addition 
to flow across the junction, which would surely favor the 
latter scenario. 



Appendix B: Transformation properties of current 
worldsheets 



Here we show that the relation between the vortex 
worldsheet and the electric current Eq. ( 29 ) , 
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has the correct properties under Lorentz transformations. 
That is, the worldsheet J^^ is a Lorentz tensor and the 
electric current is a Lorentz vector. 

Start out from a static vortex line in the z-direction. 
The only non-zero component of the vortex worldsheet 
is JTjY and the current is = (cp, 0,0,^2). Perform a 



rotation in the xz-plane over angle a. The transformed 
fields are, 
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It is easily verified that relation Eq. (Bl) holds. 

Next, from the static line in the z-direction, perform 
a Lorentz boost in the a:-direction with velocity v. The 
Lorentz transformation matrix is. 
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Here 7 
are, 



^1— u2/c2 



and P — v/c. The transformed fields 
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The worldsheet is no longer pointing 'straight up' in the 
time direction, but is tilted with a non-zero component 
J'xz, denoting the motion in the x-direction of the vor- 
tex line along z. For an Abrikosov vortex this component 
would correspond one-to-one with an electric field in the 
y-direction. But as we see, it now leads to a component of 
the electric current in the x direction. If we compare the 
two currents Eqs. (B3| and (B6), then given a 4-current 



w^, we cannot say whether it corresponds to a stationary 
or a moving vortex line. We need additional informa- 
tion, such as time derivatives. Relation Eq. 

A ^'==1 = ^iPcpwz = -yPwz = 



J', 



(Bll holds 



because — 

cp I 

Now look at a Lorentz boost of the rotated current, 
which already has a component in the x-direction. The 
transformed fields are. 
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For the boost of Jj^, , this follows from J/^ 



J[^. This component is invariant under boosts in the 
x-direction, because the contraction in the spatial x- 
direction is 'compensated' by the dilation in the temporal 
direction. The surface area of elements in the xt-plane 
is therefore unchanged. For the corresponding current 
component this follows from, 



(7U'tdf' + 7/3w^dx') A {P^w'^At' + 7w^da;') 
-i^w[w'^dt' A dx' + -I'^p'^Ax' A At' 
(7^ - 7^/32)w>;dt' A Ax' 



^w'^At' hw'^Ax' . 



(B9) 



Here we used the antisymmetry of the wedge product A. 
For the verification of the xz-component we compute, 



v'^Ax" A w'lAx" = il3-fw'tAt' + jw'^Ax') A w'^Az' 

= {jS^cpAt + 7 sin a WzAz) A cos aWzAz 

— P^cp cos aWzAt A Az 

= /37 cos aJtzAt A Az. (BIO) 



This agrees with Eq. (B7). 
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tions can never lead to a sign-free ground state, as is re- 
quired for a many-boson wave function. This was pointed 
out most recently by P.W. Anderson in Ref. US] 
The name "t ype -II Mott insulator" was used earlier by Lee 
and KivelsorpSl to denote insulators that go continuously 
into a conductive regime by doping. Here we condisder 
the type-II phenomenon of Abrikosov-like vortices in the 
undoped, bosonic case. 



